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Abstract 

An explanation is given for the observed magnetic- 
field dependence of the low-temperature specific 
heat coefficient of YI54AS3. It is based on a re- 
cently developed model for that material which can 
explain the observed heavy-fermion behaviour. Ac- 
cording to it the Yb 3+ -ions are positioned in a net of 
parallel chains with an effective spin coupling of the 
order of J w 25 K. The magnetic-field dependence 
can be understood by including a weak magnetic 
coupling J' between adjacent chains. The data re- 
quire a ratio J'/ J « 10~ 4 . In that case the exper- 
imental results can be reproduced very well by the 
theory. 



The low-carrier compound Yb4As3 which belongs 
to the Rare-Earth pnictide family has been sub- 
ject of extensive experimental investigation caused 
by its intriguing low-temperature anomalies [1, 2, 
3, 4, 5, 6, 7, 8, 9]. It shows evidence for a 'heavy- 
fermion' type behaviour below a characteristic tem- 
perature T* ~ 40 K, notably a large 7-coefficient of 
about 200 mJ/mol • K 2 for the specific heat. On the 
other hand the carrier concentration is extremely 
small, approximately 10 -3 holes per formula unit 
below T* . This excludes a conventional Kondo- 
mechanism as origin of the low-energy excitations 
leading to the large 7-value. 

The scenario proposed in Ref. [10] is quite dif- 
ferent from that of a Kondo system. The crystal 
has the anti-Th3P4-structure, and the space group 
is I43d. The Yb-ions occupy four systems of par- 
allel chains oriented along the space diagonals of 
a cube (body-centred cubic rod packing) [11]. A 
structural phase transition at room temperature 
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changes the high-temperature homogeneous mixed- 
valent state of Yb 4 As 3 (Yb 3+ : Yb 2+ = 1:3) 
to a charge ordered state where the 25 % Yb 3+ - 
holes are confined to one of the four chain systems, 
e.g, parallel (lll)-chains. The particular ordering 
of Yb 3+ -ions was confirmed clearly by perturbed 
angular-correlation (PAC) measurements [6]. This 
state has almost half-filled hole bands which are 
subject to strong on-site Coulomb correlations, ef- 
fectively described by a t-J model. Therefore the 
low-energy excitations essential for the specific heat 
are the magnetic excitations of an antiferromag- 
netic quasi-one-dimensional l S = i' Heisenberg 
chain (S is the pseudospin corresponding to the 
lowest Kramers doublet of Yb 3+ ). This model pro- 
posed in Ref. [10] was beautifully confirmed by 
the inelastic neutron-scattering results of Kohgi et 
al. [7]. They show that low-temperature magnetic 
excitations in Yb4As3 are very well described by 
the one-dimensional spin wave spectrum, uj(q) = 
(tt/2) Js'm(aq), of des Cloizeaux and Pearson [12] 
where q is the wave number in (lll)-direction and 
a the distance of Yb 3+ ions along the chain. This 
also leads naturally to a linear specific heat for 
k^T I J <C 1. From the maximum observed spin- 
wave energy of 3.8 meV ~ 40 K at q = ^(n/d) one 
obtains J — 25 K in good agreement with the en- 
ergy scale T* estimated from the value of 7. 

However there is still one major challenge to the 
validity of the model proposed in Ref. [10]. Con- 
cerning the dependence of 7 on an applied mag- 
netic field H, one would expect deviations from 
7(H = 0) of the order (H/J) 2 , i.e., about 10~ 2 at 
4 T uniformly for all temperatures with k^T / J < 1. 
Experiments by Helfrich et al. [2] revealed a quite 
different behaviour: Above 2K changes in 7 are 
indeed small. In contrast, at 0.5 K the field depen- 
dence of 7(H) is dramatic. Already at 4T its value 
is strongly suppressed in conflict with the above ar- 
gument. This indicates that an additional energy 
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scale much smaller than J (25 K) must be present 
in Yb4As3. 

It is the purpose of this letter to show that the 
spin-chain model for Yb4As 3 is well capable of ex- 
plaining the unexpected results of Hclfrich et al. 
For the present purpose we may neglect deviations 
from perfect charge ordering or half filling which is 
only important for the transport properties. The 
additional energy scale which the results in Ref. [2] 
suggest is associated with anisotropy gaps of the 
magnetic excitations caused by interchain coupling 
J' of the parallel Yb 3+ -spin chains. The complete 
excitation spectrum of a trigonal antiferromagnet- 
ically coupled chain system is calculated and it is 
shown how the linear specific-heat term and its field 
dependence arises. A comparison with the data of 
Hclfrich ct al. allows to estimate the magnitude of 
J'. Proposals for magnetic resonance experiments 
are made. 

It is well-known that in one spatial dimension, 
there cannot be any long-range spin ordering, but, 
according to calculations on the spin- 1/2 Heisen- 
berg chain [13, 14], the spin-spin correlation func- 
tion decays slowly, (S(r) • S(0)) oc (— l) r y / log(r)/r, 
i.e., the correlation length is infinite. We therefore 
model Yb4As3 as if the spins were antifcrromag- 
netically ordered and apply linear spin- wave theory, 
neglecting disorder-induced corrections. 

We transform the crystal coordinate system of 
the charge-ordered state to a new basis with a 
hexagonal unit cell where the (lll)-direction now is 
oriented along the c-axis, i.e., it transforms to the 
(OOl)-direction. Taking into account the interchain 
coupling by a parameter A, we describe Yh^Ass 
with a Hamiltonian of the form 



+ ^(d^-h.Si), 



(i) 



where A = J' / J is the ratio of the antiferromag- 
netic exchange constants perpendicular and paral- 
lel to the c-axis, and h = .9 c ffMBH/ J is the applied 
magnetic field. All energies are in units of J. We 
parameterize the anisotropic exchange of the sys- 
tem by introducing a parameter d > describing a 
single-ion easy-plane anisotropy. The symbols (ij) 
and (ij)' imply summation over nearest- neighbour 
bonds along and perpendicular to the c-axis. 

Starting point of the calculation is the classical 
Neel-state. Perpendicular to the c-direction, the 
spins are situated on a trigonal lattice-structure, 
subdivided into three interpenetrating sublattices. 



These sublattices are stacked antiferromagnetically 
along the c-direction. The ordered system thus 
contains six sublattices £(L,c), numbered with 
L = 1,2,3 within one plane perpendicular to the 
c-axis, and with c = 1,2 in c-direction. In or- 
der to describe the Neel-state properly, we per- 
form coordinate transformations on each sublattice 
£(L,c) such that each Sj points into the local z- 
direction. Parametrizing these transformations by 
polar and azimuthal angles ^l,c and @l,c, we apply 
a Holstein-Primakoff transformation on the spins 
and retain only the constant and bilinear parts of 
the Hamiltonian. After Fourier transforming H, we 
diagonalize it by a Bogoliubov-transformation. The 
diagonal form of H then reads 



H 



+ 2S E E w ^(q) (al Lc a g , Lc + M 

L,c q 



(2) 



H oc (S, h) contains constants arising from a rear- 
rangement of the bilinear terms of H, and a q x c 
( a q Lc) * s the Fourier transform of the Holstein- 
Primakoff boson dj (a\) on the sublattice C(L,c). 
i?Neei oc (S 2 , h ■ S) contains the classical part of H , 
i.e, the exact Hamiltonian for S — > oo. Minimizing 
it with respect to the angles $ l c and 6 Lc yields the 
ground-state spin configuration of (2) [15, 16]. The 
problem of quantum corrections to the ground-state 
energy may be important for some cases [17, 18]. 

If the applied field h does not break the rota- 
tional symmetry of (2) in the limit h = 0, i.e, 
for h parallel to the c-axis, the Neel-state is made 
up of two antiferromagnetically stacked triangu- 
lar lattices. Within these lattices, the spins of 
the three sublattices lie on the surface of a cone 
(\$Lc ~ $l>c\ = 27r/3, L ^ L') with the axis par- 
allel to the field (the c-axis) and an apex angle 

Qlc = 6 = arccos ^ 2 +rf+9/2A ) • ^ nc s P m -wave dis- 
persion is given by 



w| c (q) 



1-2 cos 2 9 + — cos e 



25 



+ (-l) c cosg z 



+ -A(l-3cos 2 e-2p L (q)) 

x ^1 + d - 2 cos 2 9 + ^ cos 9 
- (-l) c (l - 2 cos 2 9) cos<7 z 
|A ( J ' ms J (-) 
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- i i l--cos 2 e) PL (q) 



(3) 



Pi(q) = 



cos' 



(?-•') 



+ COS 



(f-*0 



COS 



2tt 



2tt 

3 ' ~ u "T 
According to Ref. [16], with a magnetic field per- 
pendicular to the c-axis, one obtains for a small 
single-ion anisotropy d a so-called umbrella-type 
spin configuration for h > hi w 2S\/d. In the 
umbrella-phase, the spins cant out of the c-plane — 
they lie on the surface of two cones with the axis 
parallel to the field direction, but different apex an- 
gles. Since in this case a simple closed-form solu- 
tion for ujL C (q) is not achievable, we diagonalized 
the Hamiltonian (2) numerically. 

Figure lc shows a typical dispersion unfolded in 
the chemical Brillouin zone for an interchain cou- 
pling A = 0.01, vanishing easy-plane anisotropy d, 
and vanishing magnetic field. In the directions par- 
allel to the c-axis, the scale is set by the intrachain 
coupling, which fits the neutron scattering data on 
Yb4 AS3 if we set J w 25 K, rescaling by a factor of 
2/ir according to the Bcthe-Ansatz result. In the 
directions transverse to the c-axis, the scale of the 
dispersion of the spin-wave modes is given by v^A. 
There are three modes becoming soft at the T-point, 
corresponding to the continuous degeneracy of the 
ground state with respect to rotations of all spins 
by the same angle around an arbitrary axis. 

Upon applying a magnetic field, this degeneracy 
is lifted, and the corresponding modes acquire a 
finite frequency at the T-point. For an arbitrary 
small but finite anisotropy d, we have still one Gold- 
stone mode if the field is applied parallel to the 
c-axis, corresponding to the independence of the 
ground state with respect to rotations around that 
axis, see Figure lb. For a field direction perpendic- 
ular to the c-axis, a gap opens in the entire Brillouin 
zone. For \fd < h < h s = 2 + 9A/2, the value of 
the gap is E g /J w \/3A. In Figures la and b, we 
take h — 0.2, corresponding to a magnetic field of 
about 6 T. 

The scenario outlined above allows for an at 
least qualitative understanding of the behaviour of 
the specific heat Cy in a magnetic field as fol- 
lows. For temperatures k^T j J <C V~A, the crys- 
tal behaves "three-dimensional", i.e., Cy oc T 3 . 
If k^T/J > \/A, the T-dependence becomes Un- 
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Figure 1: Spin- wave dispersion for the trigonal- 
chain structure in a magnetic field, (a) Field ap- 
plied perpendicular to the c-axis, (b) field parallel 
to the c-axis, (c) without applied field. The ratio 
A of intcrchain-to intrachain-coupling is 0.01. 
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Figure 2: Magnetic- field dependence of the 7- 
cocfficient. Bottom solid line: field applied perpen- 
dicular to the c-axis, top solid line: field parallel to 
the c-axis, both at k B T = 0.02 J w 0.5 K. Dashed 
line: field dependence at k B T = 0.2 J w 5 K. The 
open circles denote the data from Ref . [2] . The ratio 
A of interchain- to intrachain-coupling is 10~ 4 . 



ear, Cy w 7T, because then interchain-excitations 
are thermally populated, and only the quasi-one- 
dimensional intrachain-excitations contribute. In 
Figure 2 we show the field dependence of the 
7-coefficicnt at k B T/J = 0.02 (solid lines) and 
k B T/J = 0.2 (dashed line). In a magnetic field 
parallel to the chains, 7 is reduced by a factor of 
2. In this case, there is still the above-mentioned 
zero-energy mode giving the main contribution to 
7. If the field direction points perpendicular to the 
chains, 7 decreases rapidly as a function of h due 
to the opening of the gap in W£ C (q). This rapid de- 
crease compares quite well with the experiment — 
the experimental values are drawn as open circles 
in the plot [2]. The experiment was done using a 
multiple-domain crystal, so one would expect the 
data to fall between the two extremes, h parallel 
and h perpendicular to the c-axis. 

At temperatures k B T/J ^> y/A, the impact of 
h on 7 is almost completely suppressed, see the 
dotted line in Figure 2. Compared with the experi- 
mental behaviour of j(T) at sufficiently large fields, 
this provides us with an order of magnitude for the 
interchain coupling, i.e., we have A = O(10 -4 ) in 
order to describe the experiment. The strong de- 
pendence of the 7-coefficicnt on A is shown in Fig- 
ure 3. 7 decreases by almost four orders of magni- 
tude upon variation of A from zero to one, i.e., were 



Figure 3: Dependence of the 7-cocfficicnt at k B T = 
0.02 J w 0.5 K on the ratio A of interchain- to 
intrachain-coupling. Note the logarithmic scale on 
the ordinate. 



A > O(10 2 ), the characteristic heavy-fermion be- 
haviour of Yb4As3 would not be present. 

We conclude that the spin-chain model for 
Yb4As3 proposed in [10] and its extension discussed 
here is in excellent aggreement with the neutron 
scattering results of Kohgi et al. and the field de- 
pendence of the 7-coefficient found by Helfrich et 
al. In particular it explains very well why this field 
dependence sets in only at low temperatures. This 
is due to a new energy scale V J J' <C J set by the 
anisotropy gaps in the excitation spectrum which 
are caused by the interchain coupling. From the 
behaviour of 7(H) one would obtain J' / J ~ 10~ 4 . 
In this case it would however be impossible to di- 
rectly observe the branches of excitation propagat- 
ing transversely to the chain direction (Fig. 1) by 
neutron scattering. We therefore suggest that mag- 
netic resonance methods which could determine the 
T-point excitation energies would be very helpful 
in determining J' and especially d more accurately. 
Although Yb4As3 is semimetallic this should not 
be impossible since the residual resistivity (about 
1 mOhm-cm) is quite large in this compound. 
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